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Abstract. As a droplet impacts on a granular substrate, both the intruder and 
the target deform, during which the liquid may penetrate into the substrate. 
These three aspects together distinguish it from other impact phenomena in 
the literature. We perform high-speed, double-laser profilometry measurements 
and disentangle the dynamics into three aspects: the deformation of the sub¬ 
strate during the impact, the maximum spreading diameter of the droplet, and 
the penetration of the liquid into the substrate. By systematically varying the 
impact speed and the packing fraction of the substrate, (i) the substrate de¬ 
formation indicates a critical packing fraction <j>* ps 0.585; (ii) the maximum 
droplet spreading diameter is found to scale with a Weber number corrected 
by the substrate deformation; and (iii) a model about the liquid penetration is 
established and is used to explain the observed crater morphology transition. 


Liquid droplet impact on a granular layer is very common in nature, industry, 
and agriculture, ranging from raindrops falling onto desert or soil to granulation in 
the production process of many pharmaceuticals. In spite of its commonness, the 
physical mechanisms involved in the impact of a droplet on sand did not attract 
much attention until recently HM, and the underlying physics is still largely 
unexplored. In contrast, droplet impact on a solid surface or a liquid pool has 
been studied extensively [131415] . However, a granular substrate can act both solid- 
like and fluid-like [16] and many experiments have been conducted to reveal the 
response of a granular packing to a solid object impact P2H23], where the intruder 
does not deform. Droplet impact on a granular substrate adds new challenges 
to the above: First, both the intruder and the target, not merely one of the two, 
deform during impact; second, the liquid composing the droplet may penetrate into 
the substrate during the impact and may, in the end, completely merge with the 
grains. These complex interactions between the droplet intruder and the granular 
target create various crater morphologies as reported in the literature mwm 
[see Fig. [Tj for examples] . An appealing and natural question is by what mechanism 
craters are formed and how this leads to the observed rich morphological variation. 
This is the main focus of the present work. Quantitative dynamic details, e.g., on 
the deformation of droplet and substrate in response to the impact, are necessary 
to gain insight about this issue. Previous dynamic measurements only regard the 
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Figure 1. Various crater morphologies observed in the experi¬ 
ments are categorized as: doughnut, truffle, and pancake. 

droplet spreading mm or individual splashing grains [S], and a systematical study 
of the effect of the packing fraction of the substrate is still missing. In this paper, we 
perform dynamic measurements with high-speed laser prohlometry and study the 
dependence of the dynamics on impact speed and packing fraction of the granular 
substrate. Following the discussion about the results of the substrate deformation 
and the droplet deformation and the underlying physics, a quantitative model is 
established which explains the observed crater morphologies in the end. 

1. Experimental methods 

A water droplet with diameter D 0 = 2.8 mm acts as the intruder. The water 
droplet is pinched off from a needle at rest and dropped from a certain height. The 
droplet is accelerated by gravity and impacts on the horizontal surface perpendic¬ 
ularly. The impact speed, Uq, is computed from a calibrated height-speed profile. 
Rain consists of droplets with a maximum diameter 6 mm [24] falling at their 
terminal velocity. For a droplet diameter of 2.8 mm, the terminal velocity is about 
7 m/s. In this study however, the impact speed is one of the important control 
parameters and is varied from 1.35 m/s to 4.13 m/s. The target granular packing 
consists of polydisperse soda lime glass beads with diameter d g = 70 — 110 pm and 
specific density p g = 2500 kg/m 3 . While there are studies considering the effect 
of the wettability of grains w, the wettability is not varied in this work, where 
hydrophilic grains are used. Grains are put into an oven at 90 °C for at least half 
hour before experiments, after which we let them cool down to room temperature. 
Various packing fractions, </>o, are prepared by air fluidization and a tapping de¬ 
vice [Fig. [2k] • The initial packing fraction </> 0 is computed according to the surface 
height relative to the container edge before impact. During the impact, the height 
of the substrate surface, z(x, y ), is measured by high-speed laser prohlometry with 
a typical depth resolution of 0.1 mm/pixel. The impact process is captured by a 
high speed camera (Photron SA-X2) with a frame rate of 10,000 Hz, of which an 
experimental image is shown in Fig. [2jo. For dynamic analysis, the surface height 
before impact is taken as the reference z — 0. We use two laser lines for the 
prohlometry , and, as the surface deforms, the lines are dehected. We translate 
the deflection into the height function, z(r), where r is the distance to the impact 
center in the horizontal plane, i.e., r = y/(x — x c ) 2 + (y — y c ) 2 , where x c , y c are 
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(a) 




Figure 2. (a) Setup sketch. The height, z(x,y), is measured 

as a function of coordinates on the horizontal xy plane with laser 
sheets and a high speed camera, (b) An image taken during the 
impact, where the droplet is at its maximum deformation. Two 
laser lines are used for high-speed prohlometry, where the deflection 
of the laser lines indicates the deformation of the target surface, (c) 

The laser line deflection in (b) is used to reconstruct the height, z, 
as a function of the distance to the impact center, r. The dashed 
line indicates the initial surface height, z = 0 and the center depth 
is labeled Z c . (d) Time evolution z(r,t) of the crater depth (color 
coded). 

the coordinates of the impact center. Assuming axisymmetry, the impact center is 
located by determining when the deflections reconstructed from individuals laser 
lines match. Henceforth, the height function, z(r), used for further analysis is al¬ 
ways averaged over the two laser lines. An example of z(r) and its time evolution 
are illustrated in Fig. |2](c) and (d). 

2. Substrate deformation 

We use the depth of the crater center, Z c (t) as defined in Fig. [2](c) , to characterize 
the development of the crater. A typical temporal evolution of the crater depth 
is shown in Fig. [3J In the beginning the crater becomes deeper with time. At 
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Figure 3. Crater depth versus time. The maximum crater depth 
is denoted as as Z*. A straight line is fitted on the data points of 
t < (D 0 + 2Z*)/Uo■ The slope of this fit, U s , is used to indicate the 
speed of the deformation. Here, Uq = 3.19 m/s, and (fro = 0.569. 

a certain moment, the crater reaches its maximum depth Z*, after which the 
droplet contracts and transports grains mixed with it towards the crater center. 
An avalanche subsequently occurs at a longer time scale. Both of these effects 
tend to decrease Z c , i.e., make the crater shallower. Here, we focus only on the 
early stage of Z c (t ) and quantify its evolution by measuring two quantities: the 
initial speed of the deformation, U s , and the maximum crater depth, Z* [25]. The 
behavior of these two quantities is essential to understand not only the response of 
the granular substrate, but also the deformation of the droplet and the formation 
of various crater shapes. 

2.1. Deformation speed U s . The initial deformation speed U s is determined as 
the slope of a linear fit of Z c (t) within a time duration t < ti, np , where = 
(D 0 + 2Z*)/Uo denotes the impact time scale. [Fig. [3]. The ratio of the impact 
velocity and this deformation speed is plotted in Fig. [4] for all experiments. This 
ratio is largely independent of impact speed which implies that it is an inherent 
property of the granular substrate. On the other hand, the dependence on the 
packing fraction 0o indicates a transition around 0* ~ 0.585. 0 Here, we introduce 
a simple scenario to elucidate this transition. Upon impact the droplet - with 
a mass M d - accelerates a certain amount of grains in the vertical direction to 
the initial deformation speed U s , while by vertical momentum conservation the 
momentum lost by the droplet is proportional to M d (U 0 — U s ). Thus, the mass 
of grains accelerated by the droplet is (Uq — U s )/U s times larger than the mass 
of the droplet. This ratio, plotted in Fig. [4] therefore indicates the amount of 
substrate material involved in the deformation dynamics and the yield stress of 

‘Although the measurement of U s is merely the first approximation of the deformation speed, 
we confirm that (f>* is not altered by reducing the fitting region to half its size. 
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Figure 4. Main plot: The ratio of the impact velocity Uq and 
the deformation speed of the substrate U s as a function of packing 
fraction 0 q. The Y-axis, (Uq—U s )/U s , indicates the amount of grains 
involved in the deformation M s normalized by the droplet mass 
(see text), where U s is the speed of the substrate deformation [Fig. [3], 

Inset: the volume change of the substrate, AV C , for experiments with 
the same impact speed U 0 = 3.75 m/s, but three different packing 
fractions, 0 O = 0.571, 0.588, 0.602, from bottom to top. 

the substrate, i.e., the larger this ratio is, the harder the substrate is. With 
this interpretation, the observed transition at 0* in Fig. [4] is reminiscent of the 
penetration force transition which points to the dilatancy transition [26]. 

The dilatancy transition expresses a peculiar phenomena that shearing a gran¬ 
ular packing above a threshold packing fraction results in dilation (i.e., expansion 
along the perpendicular directions), whereas a loose packing is just compactified 
under such a perturbation B3- Though the word ‘dilation’ describes a change of 
the volume, previous studies have shown that the dilatancy transition is accom¬ 
panied by a force response transition |2T,i22,;26,28]. As explained above, this force 
transition feature is already encapsulated by the factor (Uq — U s )/U s represent¬ 
ing momentum transfer, and the volume change of the substrate is supposed to 
give a transition at 0* as well. We compute the volume change as the integral, 
AV C = 7T f z(r)dr 2 , for each frame. The result is plotted in the inset of Fig. [4] for 
three experiments with different 0o but the same impact speed Uq. The volume 
change is negative (i.e., the substrate compactifies) for the loosest substrate and 
indeed increases to positive (indicating dilation) with increasing 0n. Based on 
the above evidence we define a critical packing fraction 0* 0.5851] The grains 

'Note that these measurements pertain to the initial deformation of the bed, i.e., before liquid- 
grain mixing becomes important. The dilation and its transition reported here are inherent 
properties of granular packings, independent of the intruder properties. This is different to the 
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underneath the droplet intruder are forced to rearrange to </>* during the impact 
whatever the initial packing fraction </> 0 is. This value will be used to model the 
mixing between liquid and grains. 


2.2. Maximum crater depth Z*. In contrast to Fig. [4} the maximum crater 
depth, Z* plotted in Fig. [5^, depends on both the initial packing fraction, 0 O , 
and the impact speed of the droplet, Uq- It has been shown above that a dense 
substrate is more difficult to deform, therefore it should come as no surprise to see 
that the maximum crater depth Z* decreases with (ft 0 . The increase of Z* with the 
impact speed U 0 is also anticipated: a higher impact energy generates a deeper 
crater. 

Studies on solid intruder impact on a granular substrate PTlHSi, 2Dj can help to 
further understand the impact speed dependence. For a crater created by a solid 
object, two scaling laws of its depth with impact energy are suggested based on how 
the kinetic energy, Ek, is dissipated na. Assuming that plastic deformation is the 
most dissipating process, the crater volume scales with the kinetic energy, which 
yields Z* ~ E^ 3 |181l20j . On the other hand, if material ejection absorbs most 
of the impact energy, then the converted gravitational energy leads to Z* ~ Iq/ 4 . 
The volume change in Fig. 0] implies strong plastic deformations during the impact. 
However, in the case of droplet impact, not only the substrate but also the intruder 
dissipates the impact energy via deformations. To apply the above arguments, we 
therefore first need to decide upon the distribution of Ek between the droplet and 
the granular packing. 

The droplet experiences a deceleration force. This force does work on both the 
droplet and the granular target along the total displacement, |-D 0 + Z* c . This 
work transforms the impact energy Ek = M^Uq/2 into other energy forms, such 
as surface energy of the droplet and dissipation inside the substrate. Out of the 
total displacement, \Dq is the contribution of the droplet deformation, while Z* c 
is that of the substrate deformation. Such a force, Ek/(^D 0 + Z£), indicates the 
average interaction between droplet and the substrate, from which we can estimate 
the energy absorbed by the droplet Ed and that absorbed by the substrate E s 
as the work done by this interaction force to deform the droplet and substrate 
respectively: 


( 1 ) 


E,= 


Dn 


Dn + 2 Z, 


-Ek] E s — 


Do /2 + Z\ 


-Ek 


This energy distribution is estimated using the average force between the intruder 
and the target. While it has been shown that the impact force between a droplet 
and a solid surface is time dependent [T51I29] . the subsequent analysis justifies Eq. |Tj 
by hindsight. 


dilation measured from the final crater volume [3], where the residue of the liquid-grain mixture 
plays a very important role. 
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Packing fraction O 0 


FIGURE 5. (a): Maximum crater depth Z* (non-dimensionalized 

with the droplet radius D 0 ) versus packing fraction (j) 0 . The impact 
velocity, Uo, is color coded. The dependence on 0 O is stronger than 
that on Uo. (b): Maximum crater depth scaled with the cubic root 
of the energy transferred into the substrate, E s , again versus 0o- See 
text and Eq. |Tj for details. 

With measured Z* and E k , the energy distributed into the substrate E s can 
be computed for all experiments. We find that Z* oc E s is the better relation 
to collapse the dependence on U 0 [Fig. [5 Jd]. For a given (j) 0 , a power law fit gives 
Z* oc Ef with a = 0.33 ± 0.04. This scaling suggests plastic deformation as 
the main cause of energy dissipation, which stands to reason. We will show that 
the kinetic energy distribution in Eq. [1] has essential consequences to the droplet 
deformation and the crater morphology as well. 

3. Maximum droplet diameter D* d 

In contrast to solid object impact, the droplet deforms during impact [ZHUEj. 
This droplet deformation changes the contact area between liquid and grains, which 
will turn out to play an important role in determining the crater morphology. In 
this section we focus on the maximum spreading of the droplet. 

For most of our experiments, the droplet contracts after spreading, which implies 
that surface tension is the dominant stopping force of droplet deformation. In this 
regime, the Weber number We = piDoU^/a is the expected relevant dimensionless 
parameter measuring the relative importance of the kinetic energy of the impact 
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FIGURE 6. Maximum droplet deformation D* d normalized with the 
initial droplet diameter D 0 versus Weber number We = piDqUq/g 
in (a) and the effective Weber number We 1 = D 0 /(D 0 + 2 Z*)We 
in (b), where Z* is the maximum crater depth [Fig. [3] . Packing 
fraction is color coded. The dashed line in (b) indicates D d /D 0 oc 

wV 1/4 . 


droplet and the surface energy, ffere, pi is the density of the liquid and a is 
its surface tension. Previous studies about droplet impact on a solid substrate 
have proposed two models for maximum diameter, D d , that is reached by the 
droplet. One is based on energy conservation argument where the kinetic energy 
is converted to the surface energy. This model results in a scaling D d /D 0 oc We 1 / 2 
[30]. The other more recent model balances the inertial force with the surface 
tension and suggests D* d jD 0 oc We 1 / 4 [QJ. In the latter case significant amount of 
the impact energy is dissipated into internal degrees of the droplet. 

In Fig. |6jr, the maximum droplet diameter in our experiments is plotted against 
the Weber number We for different packing fractions. We observe that data with 
the same impact speed U 0 are scattered due to the variation in packing fraction, 
where a denser packing typically generates a larger droplet deformation. This can 
be understood from the dependence of the substrate deformation on the initial 
packing fraction [Fig. [2] together with Eq. CD It has been shown that a denser 
packing deforms less, thus more energy is transferred into the droplet resulting in 
a larger spreading. In the energy conservation model [30], We needs to be replaced 
with an effective Weber number WS = D 0 /(D 0 + 2 Z*) We. Furthermore, during 
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the impact the droplet experiences an average deceleration force with magnitude ~ 
E d /\D q [Eq. [Tj. Reproducing the argument in the force balance model introduced 
in ref. iH the droplet deforms to balance such a deceleration, which yields a 
scaling with WS as well. Therefore, we suggest to use We' = D ^° 2 z* ^ e f° r ^ ie 
impact on a deformable substrate, where Z* characterizes the deformation of the 
substrate. 

When the maximum droplet deformations are plotted versus effective Weber 
number in Fig. [Hfa, we find that the data collapse onto a master curve. A power 
law fit gives We^ ' , which implies that our experiments are well described 

by the force balance model of Clanet et al [T3]. Previous studies on this topic 
rather used the traditional Weber number and reported smaller scaling exponents 
no larger than 1/5 [2115]. This had been interpreted either as a viscous effect [5] or 
as a result of the density ratio of the liquid and sand [2]. However, liquids with the 
same viscosity and density but different surface tension scale differently [5j. We 
speculate that, other than viscosity, the impact energy distribution in Eq. [T]may 
be helpful to understand those inconsistent observations in the literature. The 
mechanism leading to Eq. [Tj is that the impact energy E k is distributed according 
to the relative stiffness of the intruder and the target, i.e., the balance goes always 
towards the most easy deformation. The extreme case is that either the intruder 
or the target is undeformable, in which case the energy quickly transfers into the 
deformable medium. When both intruder and target are deformable, the impact 
energy dissipates more into the ‘softer’ one. For instance, while impacting on the 
same substrate, a liquid droplet with a smaller surface tension is easier to deform, 
and in consequence the substrate behaves more solid-like, and We 1 / 4 is more likely 
to be recovered. This could explain the different scaling exponents found for liquids 
that only differ in surface tension in ref. [5] The computation of Eq. |T] requires 
the maximum crater depth, Z*, which is measured for the first time in the present 
experiments. 


4. Liquid-grains mixture and crater morphology 

Until so far we considered the deformation of the target and that of the intruder 
separately. However, they are ‘miscible’ as well, i.e., the liquid that composes the 
droplet penetrates into the granular substrate. This penetration results in a liquid- 
grain mixture which needs to be taken into account in order to understand the 
crater morphology mm- By changing impact speed Uq and the initial packing 
fraction </>o the crater morphology is found to vary systematically. According to 
the residue shape the observed craters are categorized into three groups: dough¬ 
nut, truffle, and pancake [Fig. [Tj. The residue consists of a mixture of liquid and 
grains, which is referred to as ‘liquid marble’ in ref. |9. The discrimination between 
doughnut and truffle is vague, as while changing the impact parameters one shape 
seems to be continuously transformed into the other, whereas the transition from 
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truffle to pancake is more abrupt. In this section, a model about the liquid pene¬ 
tration is established together with the knowledge already obtained in the previous 
sections. We will use this model to explain the various crater morphologies. 

The doughnut and truffle residues are formed from the grains mixed with the 
liquid which are transported towards crater center by the droplet contraction. If 
the penetration during the impact is little, after the contraction, pure liquid con¬ 
centrates at the crater center surrounded by the mixture of liquid and grains. After 
a few hundreds of milliseconds, the liquid in the center penetrates into the sub¬ 
strate due to gravity and capillary force leaving a dimple, i.e., creating a doughnut 
residue. More penetration results in a smaller amount of pure liquid at the center, 
such that the residue gradually becomes a truffle. With even more penetration 
the droplet is hardly able to contract, leaving a flat residue, the pancake shape. 
Therefore, knowing the amount of penetration of liquid between grains is crucial 
to understand the morphology phase diagram. 

The transition from doughnut to truffle is continuous, while a sharper transition 
between truffle and pancake is observed. We estimate the volume of the liquid- 
grain mixture to characterize the penetration amount and define a threshold vol¬ 
ume V* to separate doughnut/truffle and pancake regime. When this threshold 
mixture volume, V*. is reached, the droplet loses the high curvature edge that 
promotes contraction. This happens when the volume of the mixture in the sub¬ 
strate, V*, equals the volume of pure liquid above the substrate, Vi — V*(l — 0), 
where Vi is the initial droplet volume [Fig. [7| inset]: 


( 2 ) 


V* 


Vi 

2 - 0 ' 


If the time scale to reach this critical mixing volume, t m i X , is shorter than the 
impact time scale t imp = ( D 0 + 2Z*)/U 0 , the droplet contraction is suppressed, 
and pancake shapes are observed. Otherwise, the droplet is able to contract, 
mixed with grains, and forms a doughnut/truffle residue. 

To quantitatively examine the above picture, we need to formulate the mixing 
progress. Describing the granular substrate as a porous medium, we start with 
Darcy’s law: 


(3) 


Q kP 

_ _ 

A /iL 


Here, the permeability of the granular packings, k, is defined by Carman-Kozeny 
relation k = (1 — 0) 3 d 2 /(18O0 2 ), where Q is the volume flux into the porous 
substrate, P the driving pressure, /r the dynamic viscosity of the liquid, L the 
penetration depth of the liquid into the porous substrate, and A the contact area. 
In Eq. [3] we estimate the pressure gradient inside the sand as P/L. Furthermore, 
the conservation of liquid volume calls for 


Q 

A 


(1-0) 


dL 

dt 


(4) 



RAINDROP IMPACT ON SAND: A DYNAMIC EXPLANATION OF CRATER MORPHOLOGIES 


Here, we only consider L as a function of time. The total liquid volume penetrating 
into the substrate is given by (1 — <p)AL. 

Combining Eq. [3] and Eq. 0| the penetration depth L is solved as L(t) = 
\J2 Pk/ (/i(l — 4>)) t. From here, one can define the time scale for which the volume 
of the mixture reaches the critical volume, AL{t m i X ) = V *: 


(5) 


/W*\ 2 (1 - 0)/x 
h/ 2 kP 


To apply Eq. Oa few quantities need to be evaluated and explained. The contact 
area is estimated as A = ttD^ 2 /4, i.e., the contact area at the measured maximum 
spreading diameter. The packing fraction (p is evaluated by the critical packing 
fraction 0*rather than the initial packing fraction </>o as explained in Fig. [4j The 
last missing piece is the driving pressure P. There are three candidates: the 
gravitational pressure, the inertial pressure, and the capillary pressure caused by 
the hydrophilicity of the grains that tries to pull the liquid into the bed. Since 
the droplet diameter is at the magnitude of the capillary length, gravity cannot 
significantly deform the droplet. The large droplet deformation shown in Fig. [3] in¬ 
dicates that the inertial pressure is much larger than gravity. A simple experiment 
of zero impact speed, in which we determine how fast the droplet is absorbed by 
the bed, also indicates that for the used grains capillary pressure is at least two 
order of magnitude smaller than the inertial pressure. Therefore it is justified to 
consider the inertial pressure as the driving pressure. In the analysis of maximum 
droplet deformation and the introduction of the effective Weber number WS it 
was already implied that the average deceleration force experienced by the droplet 
equals E d /^D 0 [Eq. |T]. The inertial pressure is evaluated in consequence as: 


( 6 ) 


E d MiUq 

|Do A A(D 0 + 2 Z*) 


With all the quantities in Eq. [5] now defined, we can finally compare the time 
scales, t iinp and t mix , in Fig. □ 

Three features from the experimental phase diagram are recovered in Fig. [2 
(i) the transition is dominated by impact speed; (ii) a denser packing induces 
more mixture; and (iii) the expected transition discriminating doughnut/truffle 
and pancake morphology around Uq ~ 3.2 m/s is indeed indicated by the condition 
timp = tmix■ We emphasize that the inertial pressure driving penetration during 
the droplet spreading is considered to dominate the crater morphology transition, 
rather than gravity during the recession as was suggested in ref. 0} Therefore, 
the reference time scale is the characteristic period of the inertial pressure, ti mp , 
rather than the contact time scale \ZpiDq/ct j9]. Finally, it is worthy to note that 
according to Eq. [5] and Eq. [6] the ratio of the time scales t irniP /t m i X is independent 
of the substrate deformation Z*. To apply and check the model, only D* d needs to 
be measured. 
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FIGURE 7. Phase diagram of the crater morphology. On the y- 
axis we plot the ratio between two time scales, namely the impact 
time scale ti mp and the mixing time t mix . The horizontal dashed line 
indicates a ratio of 1, where the transition is expected. On the x- 
axis we have the impact speed Uq. The various crater morphologies 
are labeled with different symbols (see legend), whereas the color 
of the symbols indicate initial packing fraction. The inset shows a 
schematic of the droplet for which V = V* at maximum expansion. 


5. Conclusion 

In this paper, we study the deformation of the droplet and the granular packing 
during impact. The deformation of the granular substrate decreases with the 
initial packing fraction [Fig. [5], however, the deformation speed exhibits evidence of 
dilation and defines a critical packing fraction (ft* « 0.585 in Fig. [4] The substrate 
deformation introduces an impact energy distribution, Eq. |T] between the intruder 
and the target. An effective Weber number, We\ is defined accordingly, and the 
droplet maximum deformation, D * d , is shown to be consistent with a scaling law 

D*d ~ We* 111 in Fig. [6j which suggests a scenario where surface tension balances 
inertial pressure. Finally, based on the results of (ft* and D d , a model is established 
to describe the penetration of the liquid into the substrate. This model evaluates 
the competition between the penetration time and the impact time, from which it 
is able to explain the observed morphology transition between the doughnut/truffle 
and the pancake crater shapes [Fig. [7]. The energy distribution as given in Eq. |T] 
which is essential to understand the deformation of the target and the intruder, 
is estimated with the average interaction force. Its validation needs time resolved 
measurement of the interaction force. 

We acknowledge the discussion with X. Cheng. This project is supported by 
FOM and NWO through a VIDI Grant No. 68047512. 
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